Vitreous silica distends in helium gas: acoustic vs. static compressibilities 
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Sound velocities of vitreous silica are measured under He compression in the pressure range 0- 
6 GPa by Brillouin light scattering. It is found that the well-known anomalous maximum in the 
pressure dependence of the compressibility is suppressed by He incorporation into the silica network. 
This shows that the elastic anomaly relates to the collapse of the largest interstitial voids in the 
structure. The huge difference between the static and the acoustic compressibilities indicates that 
the amount of incorporated helium still increases at 6 GPa. 



Silica is important for countless reasons. In its vitreous 
form, 7j-Si02, it finds a considerable number of technical 
applications. Further, it is a prototypal glass-forming 
material of fundamental interest. Thus, u-Si02 has been 
extensively studied over large pressure (P) and tempera- 
ture (T) ranges [11-01 • Several anomalous physical prop- 
erties were discovered. In particular, the compressibility 
X does not show the P and T dependencies observed in 
most solids. As P is increased at room T, x first increases 
to reach a maximum around 2 GPa and it decreases there- 
after [H, Q. The variation of x with T is also anomalous. 
It decreases as T increases above ambient up to the glass 
transition [lol - [l^ . The relation of these anomalies to the 
structure of w-Si02 remains debated. 

Two different models, supported by numerical simula- 
tions, have been proposed. In a first one, the changes 
of X with T and P are related to structural modifica- 
tions associated with rotations of Si— O— Si bonds in six- 
membered rings, similar to those occurring in the phase 
transition of cristobalite [H, • This was termed a lo- 
cal progressive transition. Alternatively, in a model based 
on network flexibility, x is maximum in an intermediate 
range of P over which volume changes can be accommo- 
dated by low-energy buckling of the "floppy" tetrahedral 
network [l^. Both approaches rely on the existence of 
a sufficiently large free volume allowing strong structural 
modifications with T and P. The mass density of v- 
Si02 is indeed rather small compared to crystal quartz. 
Also, the fact that v-Si02 can be permanently densified 
by more than 20% by pressurization at 20 GPa [l^ 17 1 
indicates that the free volume is unusually large. 

The distribution of interstitial voids in glasses can be 
investigated by studying gas solubility. Results for v- 
Si02 were analyzed assuming a log-normal distribution 
of voids, with a diameter around 0.2 nm and a width of 
about 0.1 nm [iM[31j in agreement with numerical simu- 
lations [20j . Thermal motions in the glass open the voids 
to gas atoms. A solubility S '^0.0084 mol/cm'^/GPa was 
found for He in w-SiOs at P up to 0.13 GPa [21 1. 

Recent publications [22, 23| revealed that w-Si02 sub- 
mitted to high P in He atmosphere exhibits a surprisingly 



small change in volume. This was interpreted as a consid- 
erable reduction of compressibility [23j , leading to an ap- 
parent bulk modulus B = x^^ = -VdP/dV « 110 GPa 
near ambient P. This x can be called the static com- 
pressibility. Another B can be calculated from the lon- 
gitudinal (ul) and transverse (wt) sound velocities, and 



the density p, all known at ambient P, B = p{v- 



One finds an acoustic modulus B « 36.5 GPa. This large 
difference in compressibilities must relate to the open 
structure of 7J-Si02 allowing He to distend the network 
in the static limit. 

In this Letter, we present Brillouin light-scattering 
(BLS) measurements of the sound velocity for the lon- 
gitudinal and transverse acoustic modes as a function of 
P up to 6 GPa in He atmosphere [l^l- We show that 
the elastic moduli are relatively weakly affected by the 
presence of He, while the minimum in the acoustic bulk 
modulus is almost completely suppressed. 

The samples arc Suprasil F300 containing less than 
1 ppm OH supplied by Heraeus Quartzglass, Germany. 
We use platelets of 56 /xm thickness and about 100 //m 
lateral dimensions. The hydrostatic pressure P is applied 
in a Chervin-type diamond-anvil cell j25j with diamond 
culets of 800 /xm diameter. The pressure transmitting 
media are either He or a non-penetrating pressurizing 
medium (NPPM), namely a methanol-ethanol 4:1 mix- 
ture (ME 4:1). P is measured by the ruby-fluorescence 
technique (26j . 

BLS experiments are performed using the 514.5 nm 
line of a single frequency argon-ion laser. We use two 
different arrangements. The first employs a Sandercock 
3-1-3 pass system using two plane Fabry-Perot's (FP) in 



tandem [271 . The sample is placed in the platelet ge- 
ometry with an external scattering angle 0ext ~ 60 ° (see 
inset in Fig.[T]). In this case, the Brillouin frequency shift 
Au;/27r is given by Aw/wq = (2i'/c) sin ^^ext, where wq is 
the angular frequency of the incident light, v the sound 
velocity, and c the light velocity in vacuum. This config- 
uration allows transverse and longitudinal modes to be 
observed on the same spectrum, as shown in Fig. [T^. 
We also use a high resolution spectrometer (HRS) 
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FIG. 1. (Color online) (a) Brillouin spectrum at 3.4 GPa for v- 
Si02 in He, obtained with the Sandercock spectrometer. The 
platelet geometry shown in the inset is used. As usual, lines 
emitted by He (noted He in the figure) and by the diamond 
anvils (D) are observed in addition to the elastic lines (E) 
and to the transverse and longitudinal Brillouin lines of v- 
Si02 (TA and LA respectively), (b) Brillouin spectrum of 
the LA modes of t;-Si02 at 3.4 GPa measured with the high- 
resolution spectrometer in the backscattering geometry. The 
two dashed lines locate the periodically transmitted elastic 
peaks. The red line is a fit to the data. 



which consists of a plane and a confocal FP in tan- 
dem 2^ 2^. This HRS, employed in the backscattering 
configuration, allows measuring the frequency shifts of 
longitudinal modes with an accuracy of about 3 MHz. In 
this case, Auj is given by Aw/wq = (2nv/c) where n is 
the refractive index at the laser wavelength. The typical 
spectrum in Fig. [T]d shows the longitudinal Brillouin line 
at 3.4 GPa, strongly broadened, and two elastic peaks. 
Owing to the small free spectral range of the confocal 
FP, the Brillouin peak and the elastic peaks are observed 
at different interference orders. For both series of ex- 
periments, the frequency shifts and the linewidths are 
obtained by fitting the Brillouin lines to a damped har- 
monic oscillator convoluted with the instrumental profile. 
The profiles are corrected for the broadening introduced 
by the small sample size as well as for the finite 
collection aperture. 

The Brillouin frequency shift of the longitudinal 
mode was measured with both spectrometers at room 
T with P up to 6 GPa. Comparing the data 
at the same P in both geometries gives n of v- 
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FIG. 2. (Color online) Transverse and longitudinal sound 
velocities for ii-Si02 as a function of hydrostatic pressures in 
Helium (red triangles) and in ME 4:1 (blue squares). The 
accuracy of P measurements is 0.1 GPa. 
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"'?-Si02 then used to extract from the more pre- 
cise HRS data. The results are shown in Fig. [2] together 
with our previous results obtained with the NPPM . 
The very small scattering of data points should be noted. 
Comparing wl values in both media, three remarks are 
immediate: (i) the numerical values are quite similar, 
(m) however, the minimum at 2 GPa with the NPPM is 
strongly reduced in He, and (Hi) the slopes above 3 GPa 
are very similar. As for the transverse modes, they are 
silent in backscattering so that the only results for vt, 
shown in Fig. [2j are from measurements in the platelet 
geometry. The new data in the NPPM are very similar 
to these already known from the literature [H, 33|. As 
for the longitudinal mode, the numerical values for the 
two data sets shown are similar. The decrease observed 
between ambient pressure and 2 GPa under NPPM is 
reduced, but not suppressed. 

One of the most surprising results of Refs 22 , [2^ is 



the very small change in volume of u-Si02 when com- 
pressed under He. Such a strong effect of the pressuriz- 
ing medium on the static compressibility is usually ob- 
served in materials with an open porosity, such as e.g. 
zeolites Ul^. It results from pore filling by molecules 
of the pressurizing medium. For an homogeneous non- 
porous material, the relative change of volume V under 
P can be calculated from the sound velocities 13911 . 
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where Vq and po are the volume and the mass density at 
pressure Pq. In Fig. |31 we plot the V{P) measured for 
w-Si02 in a NPPM. The experimental data (squares) are 
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FIG. 3. (Color online) Relative change of volume as a function 
of P calculated using Eq. HI for ME 4:1 (blue solid line) and 
He (red dashed-dotted line), compared to experimental data 
from literature in ME 4:1 (blue squares [sj]) and in He (red 
triangles [1^). The red dashed line is a 2""^ order Birch- 
Murnaghan fit [S^ of the volume change in He. 



in excellent agreement with the values calculated from 
vi^{P) and vt{P) in ME 4:1 (solid line), as already noted 
by several authors (see e.g. Ref 0,11^). If Eq. ([J) is ap- 
plied to w-Si02 pressurized under He using the results of 
Fig. [21 the dashed-dotted curve is found. The disagree- 
ment with the experimental data (triangles) is obvious. 
Eq. dl]) fails because pV ^ poVo, which is a direct conse- 
quence of the He penetration. 

Turning now to the pressure dependence of the bulk 
modulus, the static value for ?;-Si02 in He can be de- 
rived from the volume chan ge. This is done using a 2^^^^ 
order Birch- Murnaghan fit [35[ of V{P) shown by the 
dashed line in Fig. [31 This allows calculating B plotted 
in Fig. [4^. It shows a linear increase from i? = 113 GPa 
at P = GPa to 135 GPa at 6 GPa. These values are 
surprisingly high. It is worth posing here to clarify the 
meaning of this measurement. The static B is obtained 
by comparing the sizes of the sample at different pres- 
sures. As He penetrates the network, a change in com- 
position results from a pressure change. Hence, the static 
measurement of B is not performed on a closed system. 
This B is the effective compressibility of an open system. 
It is different from the acoustic compressibility. The BLS 
experiment probes the elastic moduli of the composite 
as a closed system. Indeed, the frequency of the strain 
wave is high, « 35 GHz, so that the He concentration in 
the sample cannot change over the period of the acoustic 
wave. The interest in comparing static and acoustic com- 
pressibilities is that their difference relates to the rate of 
He inclusion in w-Si02. 

From the velocity measurements shown in Fig. [2l it is 
possible to calculate the variation with P of the elastic 
moduli, provided the corresponding values of the den- 
sity of the medium are known. In the NPPM case, this 
is easily achieved with Eq. [1] The results are shown as 
squares in Fig. [3^ and Fig. [3b- As already known, the 
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FIG. 4. (Color online) Elastic moduli calculated from the 
measured sound velocities in ME 4:1 (blue squares) and in 
He using the He content derived from Eq. Q (red triangles) : 
(a) bulk modulus, (b) shear modulus. The red dashed line in 
(a) is the bulk modulus derived from the Birch-Murnaghan 
fit shown in Fig. [3] (c) The relative He concentration tiho 
as a function of P, calculated using Eq. ([2]) (red line). The 
shadowed region is delimited by the He contents calculated 
using the other two hypotheses explained in the text. 



bulk modulus B decreases to a minimum around 2 GPa 
and increases thereafter 0, . The shear modulus G 
has a similar variation, with a much smaller amplitude. 
With He penetrating inside the network, the density to 
take into account in the equation of motion is that of a 
composite material made of w-Si02 and He atoms. It can 
be calculated if the change in volume of the sample and 
the amount of He are known. The latter can be eval- 
uated by various methods. Starting from the solubility 
S measured from ambient P to 0.13 GPa [2l|, and as- 
suming that this remains valid at higher P, the higher 
limit of the shadowed region in Fi g. [^ is obtained. A 
second estimate is proposed in Ref. |23| . It assumes that 
the volume difference between the sample under He com- 
pression and the sample compressed in NPPM is filled 
with He fluid in equilibrium. This provides the lower 
limit of the shadowed region in Fig. [H;. These estimates 
are much larger than the value 0.1 mole per mole Si02 
calculated at ambient pressure from the distribution of 
voids in a numerical simulation (20j . 

From the measured 'T-^fsi02(^) with P, we can also 
derive an estimate of the amount of He entering u-Si02. 
To this effect we apply the Clausius-Mossotti relation to 
the composite sample. 
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A^Si02 0^3102 + A^HcttHc 
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where Eq is the permittivity of free space. Here n — 
ri^f g;Q2 is the refractive index of the composite, A^si02 
and TVhg are the number densities of Si02 and He, re- 
spectively, and the a's are the polarizabilitics. For He, 
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we take anc = 2.602 x 10~^^ cm^, approximated as 
constant. -/Vsi02 is available from the observed volume 
change in He shown in Fig. [3] To estimate asi02: ^6 
assume that the volume calculated from B measured by 
BLS, the dashed-dotted line in Fig. [31 is that occupied 
by the Si02 skeleton in the composite. Its polarizability 
is taken equal to that of w-Si02 in NPPM at the same 
p, known from applying the Clausius-Mossotti relation 



to n 
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-Si02 



(P). Eq. © gives then ^hc = iVne/iVsio^, 
shown as a solid line in Fig. |4t. Remarkably, this line 
extrapolates the known solubility at low pressures [2l| . 
supporting this approach. It suggests that the He solu- 
bility remains constant over our P-range, leading to 1.3 
mole of He per mole of Si02 at 6 GPa. 

From these various evaluations of the He content, es- 
timates for p{P) are obtained which are used together 
with the measured velocities to calculate the elastic mod- 
uli. The triangles in Fig. and |3b are obtained using 
the p{P) values determined from n^lsioii-^) ^ ^ J^st ex- 
plained. The shadowed regions in the same figures cor- 
respond to the two limits for p explained above. These 
shadowed regions are very narrow, as the contribution of 
the exact amount of He to p is relatively minor. G is al- 
most constant over the range of P. From measurements 
of the variation of the elastic constants with P in crystals, 
it is known that compression quite generally produces an 
increase of the elastic constants associated with compres- 
sive modes, while those associated with shear modes can 
either increase or decrease 4l|, |42|. The conclusion is 



that in u-Si02 shear is quite insensitive to compression. 

The decrease of B with increasing P up to ^2 GPa 
in NPPM can be explained by interstitial voids that are 
progressively suppressed, and that this collapse is associ- 
ated to local structural changes making the network more 
compliant. Above 2 GPa, ?;-Si02 is sufiiciently compact 
to start behaving as a normal homogeneous solid under 
further compression. Then B increases almost linearly 
with P. A striking result of our study is that, with He as 
pressurizing medium, the anomalous decrease of B with 
increasing P practically disappears. B is nearly constant 
up to '^1 GPa and increases thereafter. This indicates 
that, even in the first stages of compression. He pene- 
trates in the largest voids and prevents their collapse, in 
line with the conclusions of earlier publications [221 123| . 
The structural rearrangements leading to the anomalous 
softening in a NPPM are no longer allowed, and the usual 
increase of B with P is then observed. At P larger than 
2 GPa, the compression of the network in presence of He 
is thus very similar to that observed with a NPPM. 

In summary, our results demonstrate that the minima 
observed in B and G when P is increased with a NPPM 
are suppressed when the compression is made under He. 
The B minimum observed in a NPPM owing to the inter- 
tetrahedral flexibility is to a large extent hindered by the 
penetration of helium atoms in the silica network. It 
would be interesting knowing from numerical simulations 



how the presence of He changes the dynamics of Si— O— Si 
flips, or of the "soft modes" , in the two models described 
in the introduction 13- Hj. Another important result is 
that the acoustic bulk modulus is less than half the static 
one. It must be stressed that the former is measured at 
constant He concentration, owing to the high frequency 
of the probed acoustic waves, while the later describes 
a composite in which the amount of He varies. In an 
open system, the static measurement underestimates the 
compressibility by a term due to the gas charging the 
composite. The large difference found here demonstrates 
that He continues penetrating w-Si02 up to our highest 
investigated P. 
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